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INTRODUCTION 


Backgro\ind 

The  U.S.  Army  Engineer  Division,  New  England  (NED)  and  the  Cold 
Regions  Research  and  Engineering  Laboratory  (CRREL)  have  been  involved 
ifi  the  Landsat  Data  Collection  System  (DCS)  and  Landsat  imagery  analysis 
since  the  laimch  of  ERTS-1,  now  known  as  Landsat-1,  in  Jvily  19T2  (Cooper 
et  al.  1975).  During  the  Landsat-l  experiment  CRREL  participated  in  the 
ICS  studies  by  developing  sensor  interfaces  for  the  Leuidsat  data  collec- 
tion platforms  (DCP's)  and  evaluating  system  performance  of  DCP  installa- 
tions. During  the  last  two  years  (1975-1977)  of  the  Landsat-2  program 
CRREL  was  involved  in  the  digital  processing  of  the  Landsat  computer 
compatible  tapes  (CCT's)  and  in  sensor  interface  development  for  the 
DCP. 


The  primary  objectives  of  the  NED/CRREL  Landsat-2  imagery  analysis 
were:  l)  to  evaluate  mapping  the  extent  of  snow,  and  2)  to  determine  if 
a relationship  exists  between  the  water  equivalent  of  the  snowpack  and 
the  radiance  obtained  from  Landsat  digital  data.  If  a general  relation- 
ship is  identified,  the  correlation  of  water  equivalent  in  the  snowpack 
with  multispectral  signatures  developed  from  the  CCT's  may  be  useful  in 
estimating  the  amoimt  of  spring  water  mmoff . Secondary  emphasis  in  the 
imagery  analysis  part  was  on  delineation  of  wetlands  and  flood  waters  in 
New  England.  Other  secondary  emphasis  in  the  total  Landsat-2  program 
was  on  developing  and  evaluating  sensor  interfaces  for  use  in  obtaining 
environmental  and  hydrologic  data  in  near-real  time  using  the  Landsat 
data  collection  relay  system. 

Symposia,  conferences  and  other  meetings  where  the  results  of  the 
Landsat-2  program  have  been  presented  are  listed  in  Appendix  A. 

Hydrologic  parameters 

A coordination  committee  comprising  personnel  from  NED,  CRREL  and 
the  University  of  Connecticut  originally  reviewed  the  hydrologic  para- 
meters that  would  affect  reservoir  operation  and  management.  The  most 
significant  hydrologic  parameters  influencing  reservoir  operation  and 
management  are: 

Snow  cover  (areal  extent,  water  equivalent  correlation) 

Soil  moisture  regimes 
Wetlands  delineation 
Slope /topography 
Ice  cover 
Flooded  areas 

The  hydrologic  parameters  selected  for  detailed  Landsat  imagery  suialysis 
using  the  Landsat  digital  data  were  snow  cover  and  delineation  of 
wetlands  and  flooded  areas. 
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INTRODUCTION  ^ 

Description  of  the  Landsat  CCT's 

The  Landsat  satellites  (Landsat-1  and  Landsat-2)  circle  the  earth 
in  a 920-km  ( 5T2-mile ) near-polar  orhit  once  every  103  minutes , each  , 

completing  approximately  ll*  orbits  per  day.  The  mxiltispectral  scanner 
(mss)  on  each  satellite  is  a line  scanning  device  which  uses  an  oscil- 
lating mirror  to  continuously  scan  perpendicular  to  the  spacecraft  (NASA 
1976).  Six  lines  are  scanned  simultaneously  in  each  of  four  spectral 

bands  for  each  mirror  sweep,  and  radiation  is  sensed  simultaneously  by  j 

an  array  of  six  detectors  in  each  of  four  spectral  bands  from  0.5  to  i 

1.1  um  (NASA  1976).  D\jring  image  data  processing  at  the  NASA  Goddard  ] 

Space  Flight  Center  a black  and  white  photographic  data  product  can  be  j 

produced  of  an  area  approximately  I85  km  (II5  miles)  on  a side  for  the  ! 

following  spectral,  regions:  MSS  band  h (0.5-0. 6 ym),  MSS  band  5 (0*6- 
0.7  ym),  MSS  b8.nd  6 (0.7-0. 8 ym)  and  MSS  band  7 (0. 8-1.1  ym).  This 
information  is  also  available  in  digital  form  on  a CCT  and  can  be 
obtained  from  NASA. 

The  standard  Landsat  CCT  was  computer-processed  to  produce  a | 

geometrically  corrected  tape  with  observations  transformed  to  a UTM  i 

(Universal  Transverse  Mercator)  projection.  This  geometrically  cor-  I 

rected  CCT  comprises  21*32  scan  lines  with  each  scan  line  covering  3200  ] 

pixels*  (Ungar  1977 ) . Differing  levels  of  radiant  energy  for  each  pixel 
within  the  scene  are  registered  on  a scale  from  0 to  127  (minimum  to 
maximiam)  for  bands  k,  5 and  6 and  0 to  63  for  band  7 (Thomas  1975)* 

Description  of  the  computer  algorithm 

The  geometric  correction  of  the  digital  data  and  the  computer 
classification  algorithms  used  in  the  analyses  were  developed  at  the 
NASA  Goddard  Institute  for  Space  Studies  (GISS)  (Ungar  1977)*  The  geo- 
metric correction  provides  for  a l:2l*,000  scale  computer  printout  which 
enables  one  to  more  accurately  locate  test  sites.  The  computer  classi- 
fication algorithms  developed  for  analysis  of  the  digital  data  allow  for 
both  components  of  the  data,  each  of  the  fo\ir  wavelength  bands  and 
associated  energy  value  for  each  pixel,  to  be  evaluated  when  classifying 
the  Landsat  data  in  various  categories.  In  addition,  atmospheric  cor- 
rections have  been  applied  to  the  Landsat  digital  data  (Ungar  1977). 

The  Landsat  MSS  observation  (pixel)  may  be  thought  of  as  a point  in 
a four-dimensional  "color"  space,  where  the  values  along  each  axis 


•Picture  element,  representing  an  area  on  the  groimd  having  dimensions 
of  61  X 76  meters. 
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represent  the  radiant  energy  received  by  the  satellite  in  one  of  the 
four  bands  (Fig.  la).  Observations  lying  in  a similar  direction  from 
the  origin  in  this  four-dimensional  color  space  are  said  to  be  similar 
in  color  regardless  of  their  total  radiant  energy.  The  distance  (length) 
of  an  observation  from  the  origin  is  a measiore  of  the  total  radiance 
associated  with  that  data  point . The  algorithm  is  primarily  designed  to 
combine  observations  that  are  similar  in  color  into  the  same  classifica- 
tion category.  , There  is  provision  for  evaluating  brightness  differences 
and  for  weighting  these  differences  in  with  the  color  discrimination 
when  constructing  the  classification  categories. 

Discrimination  based  solely  on  color  is  obtained  when  the  difference 
in  direction  between  the  color  vectors  (observations)  is  examined.  If 
the  angle  between  the  observations  is  smaller  than  some  user-defined 
ariteviorii  the  vectors  are  considered  to  be  lying  in  the  same  direction 
and,  therefore,  the  observations  are  placed  in  the  same  category. 

There  are  two  modes  in  which  this  classification  scheme  may  be 
employed:  supervised  and  unsupervised.  In  the  supervised  mode  the  user 
specifies  a signature  (the  energy  distribution  in  the  four  Landsat 
bands).  If  an  observation  lies  within  a solid  angle  smaller  than  the 
user-defined  criterion,  5 , it  is  said  to  belong  to  the  category 

represented  by  the  multispectral  signature  (Fig.  lb).  Therefore,  all 
vectors  lying  within  a cone  of  angle  5 about  the  signature^  repre- 
senting category  X belong  to  category 

In.  the  unsupervised  mode  the  color  vector  corresponding  to  the 
first  observation  is  compared  with  all  subsequent  observations.  If 
color  vector  1 is  similar  in  direction  to  color  vector  2 (i.e.  66  <^6  ), 

observation  2 is  placed  in  the  same  category  as  the  first  observation 
(Fig.  Ic).  In  a similar  fashion  observations  subsequent  to  observation 
2 are  compared  to  the  second  observation  and  so  on  right  up  to  the  last 
observation.  If  in  the  process  of  constructing  categories  a member  is 
found  which  belongs  to  a previous  category,  the  new  category  is  chained 
(or  linked)  to  the  original  classification  category,  forming  one  ,)oint 
category  (Ungar  1977).  In  effect  the  unsupervised  classification  will 
form  several  categories  based  on  a criterion  specifying  maximum  color 
difference  permissible  between  members  of  the  same  category. 

In  addition  to  discrimination  based  solely  on  color,  the  GISS  al- 
gorithm provides  the  capability  of  weighting  total  radiance  differences 
into  the  discriminant  equation  for  classification.  The  percent  difference 
in  brightness  between  two  observations  is  computed.  The  calciilated 
normalized  difference  is  then  combined  with  the  color  difference  angle 
(expressed  in  steradians)  by  performing  a weighted  average  in  the  RMS 
(root  mean  square)  sense.  This  brightness-weighted  quantity  is  now 
compared  with  the  user-defined  criterion  (6  ).  Thus,  in  the  classi- 

fication process,  a relatively  small  weight?^  of  brightness  allows  very 
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a.  A color  vector  in  a four-dimensional  space. 


Band  2 


b.  Supervised  mode.  The  user-defined  criterion,  6 , defines  cate- 

gory X about  the  specified  signature  Any  coTor  vector  that  lies 
within  this  cone  belongs  to  category  X.  This  is  illustrated  for 
three  bands;  however,  all  four  Landsat  bands  are  used  in  the  com- 
puter classification  algorithm. 

Figure  1.  The  concept  of  the  four-dimensional  "color"  space  used  in  the 
computer  classification  algorithm. 
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Band  3 


c.  Unsupervised  mode.  is  similar  in  direction  to  (66  i 

and  placed  in  category  1.  is  similar  in  direction  to  B,,  arit 
placed  in  category  2.  However  Bi^  is  also  similar  in  direction  to 
B^  (category  l).  Therefore,  category  1 is  merged  with  category  2. 

Figure  1.  The  concept  of  the  four-dimensional  "color"  space  used  in  the 
computer  classification  algorithm. 


large  brightness  differences  to  disqualify  observations  that  are  similar 
in  color  from  membership  in  the  same  category,  thereby  adding  a second 
level  of  discrimination.  Discrimination  of  the  classification  cate- 
gories based  partially  on  overall  brightness  differences  plays  an 
important  role  in  the  work  discussed  in  this  report . 

Computer  data  handling  and  analysis 

A Harris  1200  remote  Job  entry  terminal  was  utilized  at  CHREL  for 
computer  data  handling  and  analysis.  The  remote  terminal  was  used  in 
the  analysis  of  the  DCP  data  cards  (hexadecimal  format)  obtained  from 
NASA  and  the  digital  processing  of  the  Landsat  CCT's. 

Data  reduction  of  the  DCP  cards  was  accomplished  by  using  the 
Infonet  computer  system  located  in  Chicago,  Illinois.  The  computer 
programs  were  developed  at  CRREL  and  used  for  einalysis  of  the  data  from 
each  interfaced  DCP  sensor.  These  data  included  temperature,  using 
thermocouples,  and  water  equivalent  from  a snow  pillow. 

The  digital  processing  of  the  Landsat  CCT's  was  accomplished  through 
a cooperative  agreement  with  NASA  GISS.  Computer  algorithms  for  the 
analysis  of  the  digital  data  were  developed  at  GISS  (Ungar  1977)-  These 
algorithms  were  accessed  using  the  CRREL  remote  entry  terminal  to  the 
main  computer  facility  located  at  GISS  in  New  York  City. 
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Literature  review 

Manual  methods  have  been  used  to  delineate  the  areal  extent  of 
snow  and  the  mean  altitude  of  the  snowline  from  Landsat  photographic 
data  products  (Barnes  and  Bowley  197^,  Meier  19T5a,  19T5b).  However,  a 
quantitative  measure  of  the  water  equivalent  of  the  snowpack  has  not 
been  obtained  from  Landsat  photographic  data  products.  Usually  the 
areal  extent  of  snow  has  been  related  indirectly  to  subsequent  watershed 
rvinoff  occ\irring  during  the  springtime  (Meier  19T5c,  Anderson  et  al. 

197^).  Also,  the  changes  in  the  areal  extent  of  snow  cover  measured  on 
Landsat  imagery  have  been  found  to  correlate  with  changes  in  water 
equivalent  recorded  by  a snow  pillow  (Anderson  et  al.  197^)-  Another 
Landsat  manual  interpretation  method  has  used  a coded  snow  cover  classi- 
fication scheme  to  accoiint  for  vegetation  cover,  density,  aspect, 
elevation  and  slope  to  map  the  areal  extent  of  snow  (Katibah  1975)* 

A snow  mapping  experiment  comparing  the  identification  of  six  snow 
cover  types  was  accomplished  using  three  image  processing  systems — 

LARSYS  Version  3,  STANSORT-2  and  General  Electric  Image-100  (itten 
1975)*  In  addition,  other  studies  have  focused  on  digital  analyses  of 
Landsat  data  in  defining  various  snow  cover  types  (Bartolucci  et  al. 

1975,  Dallam  1975,  Luther  et  al.  1975,  AlfSldi  1976).  In  these  studies 
a quantitative  estimate  of  water  equivalent  content  associated  with  snow 
cover  types  was  not  made.  In  one  case  it  was  suggested  that  spectral 
variations  within  the  snowpack  area  could  not  be  reliably  determined 
because  of  detector  saturation  problems  (Bartolucci  et  al..  1975). 

Another  study  used  simulated  infrared  Landsat  color  composites  and 
snow  course  data  to  estimate  water  equivalent  related  to  the  snowpack 
(Sharp  1975).  Sampling  units  on  the  Landsat  image  were  mapped  to  deter- 
mine the  areal  extent  of  snow.  An  estimation  of  a snow  water  content 
index  was  calculated  using  a linear  regression  equation  relating  the 
imagery  to  ground  truth  data. 

It  has  been  stated  that  remote  sensing  of  the  snow  cover  may  have 
useful  applications  since  the  magnitude  and  wavelength  of  reflectance 
vary  with  snow  types  (Mellon  1965 ).  Also,  the  albedo  is  high  for  a 
layer  of  new  snow  and  as  the  new  snow  grains  coalesce  and  coarsen  in 
texture  the  albedo  falls  steadily  (Bergen  1975).  In  addition,  a re- 
duction in  the  spectral  reflectance  occurs  from  the  combination  of 
densification  and  increased  particle  size  associated  with  aging  (O'Brien 
and  Miinis  1975).  Therefore,  it  is  believed  that  the  Landsat  CCT's  may 
contain  information  that  can  iJe  used  to  estimate  water  equivalent  in  a 
snowpack.  If  so,  then  Landsat  digital  data  can  be  used  for  estimating 
spring  ininoff  (from  a watershed  area)  more  accurately  than  is  now  possible. 
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Approach 


Selection  of  site.  The  Dickey-Lincoln  School  Lakes  Pro.lect,  Maine 
(Fig.  2),  cvirrently  being  evaluated  by  NED  for  the  generation  of  hydro- 
electric power,  flood  control  and  recreational  purposes,  was  an  ideal 
site  for  an  analysis  of  snow  cover  utilizing  the  Landsat  CCT's.  Due  to 
the  remoteness  and  inaccessibility  of  the  area  it  is  impratical  to 
establish  an  adequate  data  collection  system  for  evaluating  the  water 
equivalent  of  the  snowpack  each  year. 


Figure  2.  Site  location  map  for  the  imagery  analysis. 

The  climate  of  this  region  is  characterized  by  short,  cool  summers 
and  long,  cold,  windy  winters.  Average  annual  temperature  is  39®F 
(3.9°C)  with  extremes  of  -Uo®F  (-Uo®C)  in  the  winter  and  99®F  {37®C)  in 
the  summer  (New  England  Division,  Corps  of  Engineers  1967).  The  average 
annual  precipitation  is  approximately  91  cm  (36  in.)  and  occurs  uniformly 
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throuf:hout  the  year  with  about  30%  in  the  form  of  snow.  The  average 
annual  snowfall  is  about  25^  cm  (100  in.)  which  occurs  during  8 months 
of  the  year  (New  England  Division,  Corps  of  Engineers  1967).  Average 
winter  snow  depth  ranges  between  51  and  102  cm  (20  and  in.),  with  the 
upper  limit  exceeding  127  cm  (50  in.).  Water  equivalent  of  the  snowpack 
reaches  a maximum  in  late  March  and  usually  exceeds  25  cm  (lO  in.).  The 
geology  and  vegetation  of  the  Dickey-Lincoln  School  Lakes  Project  area 
was  previously  mapped  and  served  as  a data  base  of  site  characteristics 
in  this  study  (McKira  and  Merry  1975,  McKim  1975). 

Ground  truth  for  the  snow  cover  analysis  consisted  of  snow  depth 
and  water  equivalent  measurements  made  in  selected  snow  courses  in  the 
upper  St.  John  River  Basin  by  the  U.S.  Geological  Survey  and  the  Allagash 
Wilderness  Waterway  Agency  for  the  winter  season  of  1972-73  (U.S.  Depart- 
ment of  Commerce  1973)  (Fig.  3,  h).  The  Allagash  B snow  co\irse  is 
located  within  a mixed  forest  near  the  confluence  of  the  Al].agash  and 
St.  John  Rivers  at  an  elevation  of  about  195  m (6^+0  ft)  msl  and  has  a 
southeast  exposxire  with  gently  sloping  terrain  (Fig.  3).  The  Beech 
Ridge  snow  course  is  located  within  a mixed  forest  near  the  Frontier- 
Churchill  Road  near  Umsaskis  Lake  at  an  elevation  of  about  396  m (1300  ft) 
msl  and  has  a western  exposixre  on  gently  sloping  terrain  (Fig.  3).  The 
Ninemile  B snow  course  is  located  within  a coniferous  forest  on  the 
floodplain  near  the  USGS  gaging  station  on  the  St.  John  River  at  an 
elevation  of  about  290  m (950  ft)  msl  and  has  a northwest  exposure 
(Fig.  3). 


Snow  course  data.  The  available  cloud- free  Landsat  CCT's  of  the 
upper  St.  John  River  Basin  were  selected  for  four  seasons  and  included 
the  following  dates:  11  February  1973  (image  ID  1203-15002),  23  July 
1973  (image  ID  1365-1^593),  26  November  1973  (image  ID  lU91-1^57I?)  and 
19  April  197^  (image  ID  1635-1^5^1).  It  was  not  possible  to  obtain  a 
cloud- free  Landsat  scene  for  February  197^,  which  would  have  been 
desirable.  Therefore,  the  11  February  1973  CCT  was  selected  for  the 
snow  cover  analysis. 

The  snow  depth  and  water  equivalent  were  estimated  from  Figure  3 
for  the  three  snow  courses  (Table  l)  for  the  date  of  11  February  1973 
(U.S.  Department  of  Commerce  1973,  Li  and  Davar  1975).  These  data 
seemed  reasonable  when  compared  to  local  climatological  data  obtained  at 
Caribou,  Maine,  located  80  km  (50  miles)  east-southeast  of  the  Dickey- 
Lincoln  area.  A major  snowstorm  occurred  on  29  January  1973  between  the 
last  snow  course  measurement  (28  Janueiry  1973)  and  the  date  of  the  CCT 
(ll  February  1973).  A total  o.f  23.^  cm  (9.2  in.)  of  snow  with  a water 
equivalent  of  1.6  cm  (0.6  in.)  was  recorded  at  Caribou,  Maine.  In 
addition,  on  8 February  1973  there  was  a minor  snowfall  of  U.l  cm  (1.6  in. 
with  a water  equivalent  of  0.3  cm  (0.1  in.).  Based  on  these  d^ta,  the 
estimated  2l*.l  cm  (9.5  in.)  was  ass\nned  to  be  reasonable. 
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Figure  3 


Location  map  of  the  Allagash  B,  Beech  Bidge  and 
Ninemile  B snow  courses. 
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Table  T.  Snow  course  data  (ll  February  1973). 


Sampling 

Snow 

Water 

Cnow  course 

liOcation 

Length 

(ft) 

points 

depth 
(in. )• 

equivalent 

iin.)» 

Allagash  B 

Ur°05'N/69°0U'W 

1000 

10 

U2 

9.6 

Beech  Ridge 

1^6°36'N/69°28'W 

- 

- 

i»l 

9.U 

Ninemile  B 

»1  in.  = 2.5I* 

1*6°1*2’N/69°U3'W 

cm 

1000 

10 

38 

9.5 

The  snow  courses  were  located  on  each  of  the  CCT's  by  generating  a 
geometrically  corrected,  l6-level  grayscale  computer  printout  of  a 320 
by  256  pixel  area  (38O.6  km^  or  ll*6.9  miles^)  at  a scale  of  1:2U,000.  Each 
observation  was  assigned  one  of  16  levels  of  gray  depending  on  its 
radiance  value  in  MSS  band  7.  The  snow  course  sites  were  located  on  the 
grayscale  printouts  using  available  topographic  maps  for  orientation. 

The  computer  test  site  containing  each  snow  course  is  hO  by  32 
pixels  for  a total  area  of  6.0  km^  (2.3  miles^).  The  snow  course  was 
located  in  the  center  of  each  computer  test  site.  The  computer  algorithm 
described  previously  was  applied  to  extract  information  concerning  the 
spectral  characteristics  of  the  snow  cover /vegetation  within  the  snow 
course  computer  test  sites. 

Results  and  discussion 


Unsupervised  classifications  were  performed  on  the  three  snow  course 
sites  for  the  11  February  1973  CCT  for  a range  of  (delmax)  values 

between  0.02  and  O.OU  with  several  brightness  weightings  (for  example,  0.1, 
0.2,  0.3)  for  initial  classification  of  the  digital  data.  Computer  runs 
which  produced  large  numbers  of  categories  were  selected  so  that  several 
signatures  could  be  extracted  for  the  pixels  contained  within  the  snow 
course  areas.  This  allowed  for  eui  evaluation  of  signature  variation 
within  each  snow  course.  The  sun  elevation  angle  (23°)  of  the  scene  was 
corrected  to  zenith  to  account  for  seasonal  variations  in  irradiance. 

The  three  snow  course  computer  classification  printouts  for  11 
February  1973  are  shown  in  Figiires  5,  6 and  7 (Merry  et  al.  1977).  The 
location  of  the  snow  course  is  shown  outlined  on  the  computer  test  site. 

The  radiance  values  associated  with  each  pixel  within  the  snow  course 
are  indicated  by  the  arrows  within  the  categorization  eumury  shown  in 
Figures  5,  6 and  7.  In  this  summary  the  two  left-hand  columns  are  the 
category  symbol  and  the  number  of  pixels  (num)  within  each  classifica- 
tion category,  respectively.  The  toZ  colimin  is  a measure  of  the  varia- 
tion among  the  signatures  within  a category,  with  the  smaller  numbers 
indicating  little  variation.  The  normalized  radiances  in  each  band 
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Figure  6.  Computer  classification  printout  of  the  Beech  Ridge  snow  course  (ll  February  1973)* 
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{Bl',  B2',  B3'  and  B4')  always  total  one  and  are  listed  for  each  cate- 
f'ory.  Also,  the  true  radiances  are  listed  for  each  band  (BJ,  B2,  B3  and 
B4);  these  values  sum  to  the  total  radiance  (mW/cm^  sr*)  listed  in  the 
extreme  right-hand  column.  The  delmax  and  the  brightness  weighting  used 
in  the  unsupervised  classification  are  shown  at  the  bottom  of  the  cate- 
gorization summary. 

The  Allagash  B snow  course  computer  classification  printout  is 
shown  in  Figure  5.  The  total  radiance  of  pixels  contained  in  the  snow 

course  area  varied  from  6.93  to  T.T^  mW/cm^  sr  (categories:  E T & P 2 

6),  which  corresponded  to  a water  equivaJ.ent  of  2U.U  cm  (9.6  in.)  obtained 
from  the  snow  course  data.  An  important  observation  was  that  the  radiance 
for  MSS  band  T was  consistently  2.99  or  3.18  mW/cm^  sr,  a difference  of 
only  one  energy  level. 

The  Beech  Ridge  snow  course  computer  classification  printout  is 
shown  in  Figure  6.  The  total  radiance  of  pixels  contained  in  this  snow 

course  varied  from  5.3^  to  6.5^  mW/cm^  sr  (categories:  * . + Q W U I 0 

G)  and  corresponded  to  a water  equivalent  of  23.9  cm  (9.*+  in.)  obtained 
from  the  snow  course  data.  The  Ninemile  B snow  course  computer  classi- 
fication printout  is  shown  in  Figure  7-  The  total  radiance  for  the 
pixels  contained  in  the  site  varied  from  5*^5  to  6.8j  mW/cm^  sr  (cate- 
gories: / * W C)  and  corresponded  to  a water  equivalent  of  2h.l  cm 
(9.5  in.)  obtained  from  the  snow  coiirse  data. 

2 

The  radiance  varied  from  5.3^  to  T.T^+  mW/cm  sr  over  the  three 
sampled  snow  course  areas,  which  corresponded  to  a water  equivalent 
value  of  approximately  2U.1  cm  (9.5  in.)  (Merry  et  al.  1977)-  The  range 
in  radiance  values  from  5*3^  to  j.jh  mW/cm^  sr  may  be  attributed  to 
variations  in  vegetative  cover,  slope  and  aspect  among  these  snow  coiirse 
areas . 

The  greatest  radiance  occurred  in  cleared  areas  such  as  fields  and 
river  channels.  As  an  example,  the  snow  cover  on  the  St.  John  River 
(Fig.  U and  6)  showed  the  highest  total  radiance,  which  ranged  from 
20.23  to  as  high  as  29.21  mW/cm^  sr.  These  high  radiance  values  occurred 
in  areas  where  there  was  a minimal  vegetative  cover.  The  important 
factors  to  be  considered  in  the  snow  cover  mapping  and  assessment  of  the 
water  equivalent  analysis  based  on  radiance  values  obtained  from  the 
Landsat  CCT's  in  the  Dickey-Lincoln  School  Lakes  area  eu:e  probably  the 
vegetation  cover,  slope,  aspect,  geomorphic  position  and,  to  a lesser 
degree,  elevation. 


^Milliwatts  per  square  centimeter  per  steradian. 


\ 


15 


! 


J A preliminary  analysis  of  the  11  February  1973  CCT  using  the  GISS 

I computer  algorithm  showed  that  the  radiance  of  the  snow  cover/vegetation 

varied  from  approximately  20  mW/cm^  sr  in  non-vegetated  areas  to  less 
i than  H mW/cm^  sr  for  densely  covered  forested  areas.  Comparison  of  the 

1 digital  data  fi-om  three  snow  courses  in  the  Dickey-Lincoln  School  Lakes 

area  with  the  radiance  value  of  the  snowpack  at  these  sites  indicated 
[ that  the  radiance  of  the  pixels  contained  in  the  snow  courses  varied 

! from  5.3^  to  T.T*+  mW/cm^  sr,  with  the  average  radiance  value  being  6.h  + 

i 0.6  mW/cm^  sr.  The  water  equivalent  of  the  snowpack  for  this  range  of 

t radiance  values  was  approximately  2U.1  cm  (9.5  in.)  of  water. 


Average  multispectral  signatures  were  extracted  for  the  three  snow 
course  sites  from  the  four-band  energy  values  for  the  dates  of  11 
February  1973  and  23  July  1973  to  evaluate  seasonal  variations  in  radi- 
ance. The  multispectral  signatures  from  these  dates  are  shown  in  Table 
TI. 
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Table  II.  Average  multispectral  signatures  (mW/cm  sr) 
for  the  three  snow  coxirses. 


Snow 

MSS 

MSS 

MSS 

MSS 

Total 

Course 

Date 

5 

6 

7 

radiance 

Allagash  B 

11 

Feb 

73 

1.7i+l 

1.189 

1.209 

3.081+ 

7.223 

23 

Jul 

73 

0.518 

0.21+7 

0.603 

1.760 

3.128 

Beech  Ridge 

11 

Feb 

73 

1.053 

0.600 

0.720 

1.927 

I+.300 

23 

Jul 

73 

0.507 

0.233 

0.561 

1.705 

3.006 

Ninemile  B 

11 

Feb 

73 

1.U28 

O.53I* 

0.975 

1+.033 

6.970 

23 

Jul 

73 

0.550 

0.282 

0.591 

1.879 

3.302 

The  area  where  the  three  snow  courses  are  located  showed  relatively 
the  same  radiance  value  (3  mW/cm^  sr)  for  the  month  of  July.  This  low 
value  would  be  anticipated  due  to  the  absence  of  snow.  The  minor  differ- 
ences observed  in  the  four  MSS  bands  can  be  attributed  to  the  difference 
in  vegetative  cover.  The  Allagash  B and  Beech  Ridge  sites  have  a mixed 
forest  cover  and  their  multispectral  signatures  are  similar.  The  Nine- 
mile  B site  is  in  a coniferous  forest  and  the  four-band  multispectral 
signatures  are  slightly  different  from  those  of  the  other  two  sites. 

Figure  8 shows  the  computer  classification  of  snow  cover/vegetation 
classes  for  the  11  February  scene  of  a selected  area  (90  km^  or  35  miles^) 
near  the  confluence  of  the  St;.  John  and  Allagash  Rivers  with  the  corre- 
lative UPGS  topographic  map.  The  multispectral  signatures  (Table  III) 
for  a supervised  classification  were  derived  from  Table  II  for  the  three 

i- 
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snow  course  sites  during  February  and  another  multispectral  signature 
was  derived  from  the  four-band  energy  values  for  the  Allagash  and  Gt. 
John  River  channels.  The  four  snow  cover /vegetation  classes  are  shown 
in  Table  TV. 

2 

Table  TIT.  Input  multispectral  signatures  (mW/cm  sr) 
for  the  four  snow  cover/vegetation  classes. 


Symbol 
on  map 

1 


Symbol 

Delmax 

Multispectral  Signat 

ures 

on  map 

MSS  1* 

MSS  5 

MSS  6 

MSS  7 

Total 

radiance 

1 

0.25 

l.Ti* 

1.19 

1.21 

3.08 

7.22 

0.3 

2 

0.25 

1.05 

0.60 

0.72 

1.93 

1*.30 

0.3 

3 

0.50 

1.1*3 

0.53 

0.98 

1*.03 

6.97 

0.3 

1* 

0.25 

5.22 

U.23 

3.1*3 

7.52 

20.1*0 

0.3 

Table  IV. 

Snow  cover /vegetation  classes  mapped 

from 

the  11 

February 

1973  CCT. 

Snow  cover/vegetation  characteristics 


Water  equivalent 
(in. )* 


Mixed  forest,  195  m (6U0  ft)  msl  elevation, 

southeast  exposure,  gently  sloping  9-6 

Mixed  forest,  396  m (1300  ft)  msl  elevation, 

western  exposure,  gently  sloping  9*^ 

Conifer'^us  forest,  290  m (950  ft)  msl  eleva- 
tion, northwest  exposure,  level  9*5 

Open  nonvegetated  areas,  lowest  elevations, 
level 


•1  in.  = 2.5U 


cm 


Patterns  of  snow  radiance  values  can  be  observed  on  the  computer 
classification  printout  (Fig.  8)  which  suggests  the  interrelationship  of 
vegetation,  slope  and  aspect.  The  1 and  2 symbols  predominate,  indicating 
mixed  forest  land  at  various  elevations.  The  3 symbol  occurs  in  isolated 
areas  such  as  hilltops  and  along  the  river  channels.  The  4 symbol 
occurs  as  expected  along  the  Allagash,  St.  John  and  part  of  the  Little 
Black  Rivers.  Also,  there  were  a number  of  unclassified  pixels  (the 
dashes)  which  can  be  seen  in  Figure  8.  This  was  as  anticipated  due  to 
other  snow  cover /vegetation  classes  that  were  not  defined  during  this 
exercise.  It  is  suggested  that  these  undefined  snow  cover/vegetation 
classes  would  be  for  areas  that  have  water  equivalent  values  either 
greater  or  less  than  2U.1  cm  (9*5  in.)  and/or  other  types  of  site 
characteristics. 
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WETLANDS  MAPPING 


Literature  review 

Landsat  MSS  imagery  has  also  heen  used  to  delineate  the  extent  of 
wetlands . The  wetland  s maps  produced  have  generally  resulted  From  a 
visual  interpretation  of  the  Landsat  photographic  data  products. 

Mapping  accuracies  between  TO  and  85^  have  been  achieved  using  Landsat 
MSS  photography  (Seevers  et  al.  19T5,  Anderson  et  al.  1973a,  1973h, 

Niger  et  al.  1975,  Rehder  and  Quattrochi  1976). 

Accurate  inventories  of  wetlands  larger  than  ten  acres  were  made  in 
Nebraska  for  four  categories:  open  water,  subirrigated  meadows,  marshes, 
and  seasonally  flooded  basins.  The  inventories  were  made  by  using 
imagery  from  two  seasons  and  an  electronic  image-enhancing  system. 

Positive  print  enlargements  of  MSS  bands  5 and  7 at  a scale  of  1:250,000 
(acquired  in  the  spring)  as  well  as  band  7 (acquired  in  late  summer) 
wePe  used  to  delineate  all  wetlands  (Seevers  et  al.  1975).  Electronic 
enhancement  of  MSS  band  6 (acquired  in  the  fall)  was  used  as  an  aid  to 
differentiate  marshes. 

A wetlands  map  of  Wisconsin  was  prepared  at  a scale  of  1:500,000 
using  Landsat  MSS  bands  5 and  7 and  an  additive  color  viewer  as  a data 
ejialysis  system  (Frazier  et  al.  1975).  Wetland  areas  in  this  investiga- 
tion were  defined  as  those  which  had  enough  water  during  June  to  adversely 
affect  the  infrared  reflectance  of  growing  plants.  These  included  areas 
with  wetland  cover  types  (marsh,  sedge  meadow,  shrub-carr  fJid  lowland 
forest)  and  poorly  drained  agricultural  cropland  areas.  Tlie  primary 
criteria  for  delineation  of  wetlands  were  the  reduced  infrared  reflec- 
tance of  broad-leafed  plants  growing  in  wet  areas , the  dark  red  tone  of 
spruce  bogs,  and  the  black  color  of  organic  soil  areas  observed  when 
using  the  additive  color  viewer  (Frazier  et  al.  1975). 

Significant  changes  occurred  in  the  size  of  wetlands  because  of 
seasonal  fluctuations  in  vegetation  characteristics  and  precipitation 
(Rehder  and  Quattrochi  1976).  The  dynamic  characteristics  of  the  wet- 
lands were  not  attributed  exclusively  to  seasonal  factors,  as  signifi- 
cant changes  in  wetland  morphology  were  found  to  occirr  within  the  MSS 
bands  for  the  same  date. 

The  following  features  were  determined  from  Landsat  MSS  band  5 and 
7 imagery  enlarged  to  a scale  of  1:250,000  for  test  sites  located  in 
Maryland  and  Georgia:  a)  upper  wetlands  boundary;  b)  drainage  patterns 
within  the  wetlands;  c)  plant  communities  such  as  Spavtina  altermi flora, 
Spartina  patens,  Junaus  roemerianus;  d)  drainage  dioches  associated  with 
agriculture;  and  e)  lagoons  for  waterside  home  development  (Anderson  et 
al.  1973a,  1973b).  Mapping  at  a scale  of  1:250,000  was  adequate  for  the 
general  delineation  of  large  marshes  and  for  rather  gross  wetland  species 
associations . 


/ 
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In  addition,  digital  processing  of  Landsat  MSS  imagery  has  heen 
used  to  map  wetlands  (Anderson  et  al.  1973h,  Cartmill  1973,  Flores  et 
al.  1973,  Klemas  et  al.  1973).  Mapping  accuracies  ranging  generally 
from  78  to  99^  have  been  achieved  using  digital  processing  techniques. 

Seven  categories  of  marsh  vegetation  and  marsh  features  were 
Identified  at  an  approximate  scale  of  1:20,000  (Anderson  et  al.  1973b). 
These  categories  included  water,  sandy  mudflat,  sparsely  vegetated 
organic  mudflat,  spoil,  Iva  fvutescens,  Spartina  patens  and  Spavtina 
altemi  flora. 

Seventeen  to  thirty  spectrally  homogeneous  land  use  classes  were 
defined  in  the  Texas  coastal  zone  using  two  clustering  algorithms  avail- 
able from  the  MSA  Johnson  Space  Center  (Flores  et  al.  1973).  Many 
classes  were  identified  as  being  homogeneous  featiires  such  as  water 
masses,  salt  marsh,  beaches,  pine,  hardwoods,  and  exposed  soil  or  con- 
struction materials,  with  most  classes  identified  as  mixtures. 

Eight  vegetation  and  land  use  discrimination  classes  were  selected 
to  map  and  inventory  the  significant  ecological  communities  in  the 
coastal  zone  of  Delawai'e  (Klemas  et  al.  1973).  These  classes  were: 
Phragmites  aormtunis  (giant  reed  grass),  Spartina  altemiflora  (salt 
marsh  cord  grass),  Spartina  patens  (salt  marsh  hay),  shallow  water  and 
exposed  mud,  deep  water  (>2  m or  76  ft),  forest,  agriculture,  and  ex- 
posed sand  and  concrete.  The  Spartina  altemiflora  was  discriminated 
with  an  accuracy  of  9^-100^,  the  Phragmites  communis  showed  a classifi- 
cation accuracy  of  Q3%,  but  the  discrimination  of  Spartina  patens  was 
only  52%. 

Approach 

2 

The  site  selected  for  the  wetlands  mapping  analysis  was  a 12U-lan 
(U8-mile^)  area  of  the  Merrimack  River  estuary  (Fig.  2).  This  area  con- 
tained the  largest  variety  of  land  use  and  vegetation  classification 
units  to  be  found  in  the  Merrimack  River  Basin.  In  addition,  the 
Merrimack  River  Basin  had  been  a primary  test  site  for  the  NED-CRREL 
Skylab  Earth  Resources  Experiment  Package  (EREP)  pro.lect  (McKim  et  al. 
1975b).  Land  use  maps  prepared  from  satellite  and  aircraft  photographic 
data  products  were  available  for  this  site. 

Results  and  discussion 

The  CCT's  were  obtained  of  the  Merrimack  River  estuary  for  6 July 
1976  (image  ID  5****^-l^+082) . A grayscale  computer  printout  of  the 
Merrimack  River  estuary  was  obtained  for  the  purpose  of  locating  potential 
training  sites  for  wetlands.  In  addition,  a 1:2U,000  NASA  RB-57/RC-8 
photograph  and  a USGS  topographic  map  of  the  study  area  were  available 
for  ground  truth  comparison  to  the  computer  printout.  An  overlay  was 
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prepared  from  the  photograph  showing  the  delineation  of  the  water  and 
the  extent  of  wetlands.  Two  training  sites  were  located  on  the  north 
side  of  the  estuary  for  use  in  the  wetlands  analysis  (Fig.  9). 


70*32' 30"  70*  50' 00“  70*47*30* 


Figure  9.  Location  map  of  the  two  training  sites  used  in  the  wetlands 
mapping  einalysis. 
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A "groiind  truth"  computer  algorithm  developed  at  GISS  was  used  to 
determine  an  average  multispectral  signature  for  the  wetlands  and  the 
water  classification  categories.  This  computer  program  allows  one  to 
"tag"  certain  pixels  of  a category  into  a 32  by  ItO  array;  these  specified 
pixels  in  the  array  are  then  used  in  the  computation  of  an  average 
multispectral  signature. 

The  wetlands  overlay  that  was  prepared  from  a photointerpretation 
analysis  of  the  1:2U,000  photograph  was  placed  over  the  grayscale 
computer  printout  (scale  1:2U,000)  of  the  Merrimack  River  estuary.  The 
pixels  that  were  within  the  boundaries  of  the  delineated  wetleuids  euid 
the  water  classification  categories  were  tagged  as  the  pixels  to  be  used 
in  the  specified  array  for  the  groiind  truth  computer  program.  The 
ground  truth  computer  program  determined  the  average  multispectral 
signatures  for  the  two  specified  categories  (Table  V). 

2 

Table  V.  Average  multispectral  signatures  (mW/cm  sr) 
for  the  wetlands  and  water  categories . 


Category 

MSS  It 

MSS  5 

MSS  6 

MSS  7 

Wetlands 

0.581+ 

0.326 

0.581 

1.510 

Water 

0.551 

0.237 

O.lltO 

0.131 

The  miilti spectral  training  signatures  (Table  V)  were  used  in  a 
supervised  classification  of  the  northern  portion  of  the  Merrimack  River 
estuary.  Various  values  of  delmax  were  used  until  an  optimum  computer 
classification  map  was  obtained  of  the  wetland  areas.  The  computer 
classification  map  is  shown  in  Figiire  10.  The  symbol  1 indicates  a 
wetlands  unit,  the  symbol  2 indicates  a water  unit,  and  the  dashed  lines 
(-)  indicate  unclassified  pixels.  The  wetlands /water  photointerpreta- 
tion map  overlay  was  used  in  comparing  the  accuracy  of  the  computer  map 
of  wetlands.  A classification  acc\iracy  of  75!^  was  obtained  for  the 
wetlands  unit,  talking  into  account  the  misclassified  and  the  unclassified 
wetlands  pixels. 

The  reason  for  the  75/^  acc\iracy  may  be  the  variability  of  wetland 
species,  because  it  was  assumed  that  all  the  wetlands  contained  the  same 
vegetation  type.  If  there  was  a large  variability  in  species,  there 
would  be  different  multispectral  responses.  Also,  changes  in  moisture 
conditions  and  tidal  fluctuations  wo\ild  contribute  to  mxilti spectral 
variations.  This  would  prevent  the  wetland  areas  from  being  classified 
in  one  broad  category. 
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Input  lignMura  abl« 

Syw  Nuiti  OdmM  81  B2  B3  B4  Alb«ao  VKJ 

1 1172  0.1S0  0 58  0 33  0 58  1.51  3.00  0 20 

2 2586  0.200  0.55  0 24  0.14  0.13  1 06  0.20 


Actual  tiytatura  tatila 


81  82  B3  B4  Alt-do 

0.74  0 42  0.88  1.77  3.82 

0.80  0.25  0.15  0.12  1.12 

No.  poinn  clwifiod  ■ 3737  No.  points  undtsufisd  ■ 1383 

Figure  10-  Computer  classification  printout  of  wetlands  for  the 
Merrimack  River  estuary  (image  ID  5UUl<-lUo82) . 
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MAPPING  OF  FLOODED  AREAS 

Literature  revlev 

Landsat  MSS  data  have  been  used  for  flood  observations  because  of 
the  relatively  high  resolution,  cartographic  fidelity  and  the  near 
infrared  sensors  (Bango  1975)-  Flood  area  measurements  for  areas  of  100 
km^  (39  miles^)  or  more  have  been  made  with  less  than  a error  (Bango 
and  Salomonson  197*+). 

Flood-prone  areas  have  been  shown  to  have  multispectral  signatures 
indicating  categories  of  natural  vegetation,  soil  characteristics  and 
cultiiral  features  which  are  different  from  the  signatures  of  surroxinding 
non- flood-prone  areas  (Bango  and  Anderson  197*+).  These  differences  have 
developed  over  a period  of  time  in  response  to  increased  flooding  fre- 
quency which  enabled  the  signatures  to  be  distinguished  from  the  non- 
flood-prone  areas  (Bango  1975).  The  areas  subject  to  flooding  along  the 
Mississippi  River  were  identified  by  observation  of  various  floodplain 
indicators  such  as  natural  and  artificial  levee  systems , soil  differ- 
ences, agricultural  pattern  and  vegetation  differences,  upland  bounda- 
ries, backwater  areas  and  special  flood  alleviation  measures  in  urban 
areas. 

Landsat  imagery  has  also  been  used  to  trace  the  details  of  inunda- 
tion and  drainage  of  flood  areas  Euid  deltaic  lowlands  (Burgy  1973).  An 
overall  lightening  of  tone  observed  on  MSS  band  7 imagery  of  the  Andrus 
Island  flood  area  in  California  was  attributed  to  an  increase  of  bottom 
reflection  with  the  lowering  of  the  water  level. 

Flood  inundation  mapping  has  been  accomplished  on  MSS  band  7 
imagery  enlarged  to  a scale  of  1:250,000  acquired  one  week  to  10  days 
after  a flood  (Heillberg  et  al.  1973,  Bango  and  Sadomonson  197*+,  Morrison 
and  Cooley  1973,  Schwertz  et  al.  1976).  The  inundated  areas  showed 
shsurply  reduced  infrared  reflectance  on  MSS  band  7 because  of  standing 
surface  water,  excessive  soil  moistiire  and  stressed  vegetation  (Hadlberg 
et  al.  1973).  These  data  compared  favorably  with  flood  extent  mapping 
accomplished  on  low  altitude  aircraft  photography.  Also,  areas  affected 
by  severe  sand  and  gravel  erosion  and  sediment  deposition  were  detected 
on  MSS  band  5 (Morrison  and  Cooley  1973). 

Color  enhancement  techniques  were  used  to  produce  a variety  of 
multispectral  color  composites  at  a scale  of  1:250,000  of  flooding  along 
the  Mississippi  River  in  the  spring  of  1973  (Deutsch  et  al.  1973,  Deutsch 
and  Buggies  197*+).  Two  color  composites  of  MSS  bands  6 and  7 were 
enlarged  and  registered  to  1:250,000  scale  topographic  maps  and  used  as 
the  data  base  for  preparation  of  flood  image  maps.  Specially  filtered 
3-color  composites  of  MSS  bands  5,  6 euid  7 and  U,  5 and  7 were  used  to 
aid  in  the  data  interpretation.  In  addition,  two-color  temporal  com- 
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posites  of  pre-flood  emd  post-flood  MSS  band  7 images  were  used  in 
interpretation.  These  indicated  that  flooding  caused  changes  in  surface 
reflectance  characteristics,  making  it  possible  to  map  the  flooded  areas 
after  the  flood  waters  had  receded.  Also,  Landsat  MSS  data  have  been 
digitally  processed  to  produce  water  distribution  maps  emd  map  overlays 
that  show  the  areal  extent  of  flooding  during  the  1973  flood  of  the 

lower  Mississippi  River  (Williamson  197*+,  1975).  j 

‘ Landsat  MSS  band  7 digital  data  were  superimposed  on  a digital 
image  display  and  manipulation  system  (IDAMS)  developed  at  Goddard  Space 
Flight  Center  (Rango  and  Anderson  197*+ ) . This  enabled  a quantitative 
change  detection  analysis  for  determining  normal  surface  water  area  and 
areas  susceptible  to  flooding.  Also,  the  General  Electric  Mxaltispectral 
Information  Extraction  System  (GEMS)  was  used  to  classify  and  measure 
water  areas  according  to  differences  in  reflecteince  resulting  from 
physical  differences  in  depth  and/or  sediment  load  (Rango  and  Anderson 
197*+).  Therefore,  in  genersil,  preliminary  digital  Landsat  flood  and 
flood- prone  area  maps  have  been  produced  at  a scale  of  1:62,500;  how- 
ever, most  mapping  has  been  done  on  a regional  basis  at  a scale  of 
1:250,000. 

Approach 

The  site  selected  for  the  flooded  areas  mapping  analysis  was  the 
Franklin  Falls  reservoir.  New  Hampshire  (Fig.  2).  During  the  last  four 
days  of  Jtine  1973,  a strong,  moist  tropical  airflow  in  conjunction  with 
a stationary  frontal  system  resiolted  in  moderate  to  heavy  rain  over  much 
of  New  England.  For  example,  the  total  rainfall  was  13.0  cm  (5-1  in.) 
for  the  three-day  period  ending  at  0800  hours  1 July  1973  in  the  Frank- 
lin Falls  reservoir.  In  the  northern  portions  of  the  Merrimack  River 
Basin  this  storm  caused  the  largest  summer  flood  on  record. 

Sixty-six  percent  of  the  storage  capacity  at  Franklin  Falls  was 
utilized  in  controlling  the  flood  waters.  Large  areas  were  inundated 

for  periods  of  one  to  two  weeks  (McKim  et  al.  1975a).  This  flood  was  i 

unusual  because  of  its  magnitude,  the  extremely  high  concentration  of 
suspended  sediment  in  the  flood  waters , and  the  fact  that  it  occurred  at 

the  height  of  the  growing  season.  On  6 July  1973  a Landsat  pass  occurred  i 

over  the  New  Hampshire  area  at  peak  flood  conditions . Due  to  partial 

cloud  cover  the  entire  s\jrface  area  of  the  flood  waters  could  hot  be 

accurately  delineated;  however,  the  satellite  imagery  did  provide  a look 

at  peak  flood  conditions  in  the  Franklin  Falls  reservoir. 

Resvilts  and  discussion 

The  areal  extent  of  water  was  best  displayed  in  the  near  infrared 
band,  MSS  7-  Therefore,  MSS  band  7 grayscale  printouts  at  a scale  of 
1:2*+, 000  were  obtained  of  the  Frsuiklin  FsJ-ls  reservoir  area  for  the 


dates  of  27  October  1972  ( image  ID  IO96-I5065 ) at  low-water  stage  and  6 
July  1973  (image  ID  13U8-I506I+ ) during  flood  stage.  The  MSS  band  7 
grayscale  overprint  symbols  representing  energy  intensity  levels  were 
used  in  differentiating  the  reservoir  water  from  the  land. 

Figure  11  shows  the  Franklin  Falls  reservoir  area  on  27  October 
1972.  The  extent  of  water  is  delineated  on  the  grayscale  printout  by  a 
solid  line.  The  dotted  line  indicates  the  outline  of  the  Pemigewasset 
River  (elevation  ranging  from  98  to  110  m (320  to  360  ft))  obtained  from 
uses  topographic  maps.  Both  outlines  show  very  good  agreement  of  the 
areal  extent  of  water  at  the  low-water  reservoir  stage. 

Figure  12  shows  the  Franklin  Falls  area  on  6 July  1973.  The  extent 
of  wateo'  is  delineated  by  a solid  line  on  the  grayscale  printout  and  the 
dashed  line  indicates  the  extent  of  water  in  the  northern  part  of  the 
reservoir  area,  which  had  to  be  estimated  because  of  cloud  cover.  The 
dotted  line  shows  the  maximum  inundation  level  of  112.6  m (369.5  ft) 
for  6 July  1973  delineated  from  USGS  topographic  maps.  Again,  the 
agreement  of  the  areal  extent  of  water  between  the  computer  printout  and 
the  ground  truth  data  is  extremely  good. 

The  number  of  acres  of  water  was  quantified  for  the  27  October  and 
6 July  Landsat  scenes.  Table  VI  shows  the  total  number  of  pixels  and 
acreage  determined  to  be  water  within  the  reservoir  area  for  both  dates. 
It  shows  approximately  60%  more  water  in  the  Franklin  Falls  reservoir  on 
6 July  1973  than  on  27  October  1972,  which  was  as  expected  (McKim  et  al. 
1975a). 

Table  VI.  Areal  extent  of  water  within  the  Franklin  Falls 

reservoir  area  on  27  October  1972  and  6 July  1973* 

Pixels  Area 

Date  (no. ) (acres )* 

27  October  1972  790  1011.2 

6 July  1973  1333t  1706.2 

2 

■*1  acre  = U0U6  m 
^Estimated 
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I PART  II.  DCP  SENSOR  INTERFACE  DEVELOPMENT 

INTRODUCTION 

! Secondary  emphasis  in  the  Landsat-2  program  was  on  developing  and 

1 evaluating  sensor  interfaces  for  use  in  obtaining  environmental  and 

hydrologic  data  in  near-real  time  using  the  Landsat  data  collection  relay 
system.  The  following  sensor  interfaces  were  developed  and  tested  under 
field  conditions:  snow  pillow,  water  qu£ility  monitor,  thermocouples  and 
a tensiometer /transducer.  These  fovir  DCP  sensor  Interface  systems  are 
described  in  the  following  sections. 

i 

Snow  pillows 

i 

! During  the  summer  of  1975  two  snow  pillow  transducer  systems  were 

interfaced  to  a General  Electric  DCP  and  installed  at  Ninemile  on  the 
St.  John  River  and  at  Michaud  Farms  on  the  Allagash  River.  A circuit 
diagram  of  the  interface  system  is  shown  in  Figure  13. 

The  computer  program  used  for  decoding  the  data  is  shown  in  Appendix 
B.  A graph  of  the  water  equivalent  data  from  these  two  sites  during  the 
1975-76  winter  season  is  showr,.  in  Figiire  lU.  The  sudden  increase  in 
water  equivalent  around  early  April  for  the  Ninemile  site  cannot  be 
explained  and  is  probably  not  real. 

The  transducers  used  in  the  Interface  were  tested  under  controlled 
temperatiore  and  pressure  conditions  during  the  summer  of  1976.  The 
resulting  temperature  calibration  curve  for  the  transducers  indicated 
that  the  system  became  erratic  below  0°C  (32°F).  Therefore,  a CRREL  in- 
house  study  on  the  reliability  of  a ntunber  of  transducers  was  initiated 
and  a different  transducer  was  used  to  replace  the  original  transducer 
in  the  snow  pillow  interface. 

The  snow  pillow  transducer  system  using  the  second  model  was  used 
during  the  1976-77  winter  season.  A circuit  diagram  of  this  Interface 
is  shown  in  Figiire  15.  Two  snow  pillow  interface  installations  were 
emplaced,  one  at  the  Allagash  Falls  location  in  northern  Maine  and  one 
at  NED,  Waltham,  Massachusetts.  However,  incorrect  data  were  obtained 
from  these  two  systems  d\Lrlng  the  1976-77  season.  The  problems  encountered 
were  inadvertent  breakage  of  the  transducers  during  handling  and  shipping, 
and  unexplained,  questionable  data  telemetered  by  the  transmitter  within 
the  DCP. 

Water  quality  monitor 

A DCP  was  installed  on  the  St.  John  River  in  northern  Maine  at  the 
Dickey  Bridge,  approximately  one  mile  upstream  of  the  confluence  of  the 
St.  John  and  Allagash  Rivers.  A Martek  water  qu*d.ity  monitor  interfaced 
to  the  DCP  transmitted  the  following  water  qtiality  information:  pH, 
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data  Oita  laad  to  pin  i.  J3. 


a.  iNTERCONNECTION  DIAGRAM  — Snow  pillow  tranaducar  to  intarfaca  to  DCP 


b.  CIRCUIT  DIAGRAM  . Intarfaca 


Figure  13.  Circuit  dieigram  of  the  snow  pillow  Interface  used  in  the 
1975-76  winter  season. 
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Figure  lU.  Water  equiv41ent  of  snow  data  for  the  1975-76  winter  season. 
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Figure  15.  Circuit  diagram  of  snow  pillow  interface  used  in  the 
1976-77  winter  season. 


TtmptrQtiira,  *C  Conductivity, 


dissolved  oxygen,  river  staige,  water  temperatxire  and  conductivity.  The 
sensor  interface  was  developed  and  the  computer  program  used  to  decode 
the  data  written  during  the  Landsat-l  experiment  (Cooper  et  al.  1975* 
McKim  et  al.  1975c). 

The  water  quality  data  from  this  exercise  are  shown  in  Figure  l6. 
The  dissolved  oxygen  probe  operated  intermittently  and  the  river  stage 
measured  less  than  6l  cm  (2  ft);  therefore,  these  data  were  not  included 
in  Figure  l6.  The  data  from  the  pH  probe  between  10-17  August  probably 
did  not  fluctuate  as  indicated  on  the  graph.  The  water  qimlity  data 
compared  favorably  with  on-site  analysis  of  these  parameters  during  the 
first  week  of  operation.  The  water  quality  information  will  serve  as 
part  of  the  baseline  data  for  the  upper  St.  John  River. 
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Taermocouple  interface 


I 


A site  for  testing  a thermocouple  interface  to  monitor  air  and 
ground  temperatures  was  located  at  Sugarloaf  Mountain,  Mains  (U5°02'56"N, 
T0°23'2l"W).  This  was  the  first  time  the  thermocouple  interface  was 
test(id  under  field  conditions.  The  data  are  presently  being  evaluated 
for  accuracy.  The  emplacement  and  interfacing  techniques  developed 
d\iring  this  field  experiment  will  be  used  for  installation  of  thermo- 
couples in  Alaska.  When  validated,  this  temperatvire  measurement  system 
could  be  used  in  reservoirs  to  monitor  water  temperature  on  a daily 
basis. 

There  is  evidence  of  permafrost  on  the  upper  305  m (lOOO  ft)  of 
Sugarloaf  Mountain.  The  s\jmmit  of  the  motintain  is  veneered  with  active, 
turf-banked  terraces  which  have  moved  downslope  during  the  past  five 
years  at  a rate  of  approximately  5-1  cm/yr  (2  in./yr)  (Boms  1975). 

When  a gondola  ski  lift  was  installed  in  1967,  the  subsurface  was  found 
to  be  frozen  in  August  and  the  water  in  the  drill  holes  had  to  be  melted 
before  the  tower  foundations  could  be  installed.  In  addition,  the  two 
towers  at  the  top  of  the  momtain  moved  between  1967  and  1976  and  had  to 
be  realigned.  The  deep  pit  at  the  summit  that  houses  the  counterbalance 
weight  used  in  tensing  the  gondola  lift  line  contains  water  that  is 
frozen  year-ro\md.  The  objectives  of  this  study  were  to  evaluate  the 
reliability  of  the  thermocouple  interface  diiring  extreme  winter  weather 
conditions  and  to  monitor  the  gro\ind  temperature  to  a depth  of  30  m (lOO  ft) 
and  determine  if  permafrost  conditions  exist  at  this  latitude/elevation. 

A complete  circuit  description  of  the  thermocouple  interface  is 
given  in  Appendix  C.  The  thermocouple  interface  accepts  inputs  from  112 
copper-constantan  thermocouples.  The  inputs  are  arranged  in  l6  "banks" 
of  seven  inputs.  One  bank  of  inputs  is  recorded  with  each  update  of  the 
DCP.  A total  of  28  thermocouples  are  being  monitored  at  Siagarloaf 
Mountain  with  seven  readings  being  transmitted  in  one  Landsat  message 
(Table  VII). 

Table  VII.  Arrangement  of  the  28  thermocouples  in  four  banks. 


Thermocouple 

Depth 

Bank 

No. 

(ft)* 

I 

1,5,9,13,17,21,25 

0,8,21*, 1*0, 56, 72, 88 

II 

2, 6,10,11*, 18,22, 26 

0,12, 28,1*1+, 60,76, 92 

III 

3,7,11,15,19,23,27 

0,16, 32, 1*8,61+, 80, 96 

IV 

1*, 8, 12, 16,20,21*, 28 

1*, 20, 36, 52, 68, 81*, 100 

*1  ft  = 0.3  m 

The  temperature  measurement  range  of  the  thermocouple  interface 
unit  is  -3U  to  +32°C  (-29  to  +90 °F).  The  resolution  of  a DCP  data  word 
is  ;+  0.25°C  (^0.5°F),  or  10  pV,  and  the  copper-constantan  wire  is 
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guaxcjiteed  to  be  accixrate  within  0.5  to  0.75°C  (0.9  to  1.35°F).  There- 
fore, the  accuracy  of  the  temperature  data  is  + 0.5°C  (0.9®F).  The 
results  from  this  experiment  will  be  reported  after  the  field  inter- 
facing techniques  have  been  successfully  tested. 

Tensiometer /transducer  system  interface 

A tensiometer/transducer  feystem  has  been  successfully  interfaced  to 
the  Landsat  data  collection  system  (McKim  et  al.  1975a).  The  interface 
system  enables  moisture  tension  and  soil  volvnnetric  moisture  content 
data  to  be  obtained  in  near-real  time. 

The  Instrument  is  currently  being  tested  londer  simulated  field 
conditions  at  CRREL.  Preliminary  results  indicate  that  tension  as  low 
as  10  cm  (U  in.)  and  as  high  as  900  cm  (35^  in.)  of  water  can  be 
obtained. 

' -3 

Typical  data  for  a soil  that  has  a bvilk  density  of  1.37  g/cm  with 
a specific  gravity  of  2.63  and  a volume  of  voids  of  U7.956  are  shown  in 
Table  VIII.  In  the  initial  series  of  tests,  tensiometer  values  could  be 
obtained  that  ranged  from  about  10  to  150  cm  (1*  to  59  in.)  of  water. 
Previously  it  has  been  very  difficult  to  get  accurate  and  reliable 
numbers  on  field  moisture  tension  less  than  300  cm  (II8  in.).  It  is 
suggested  that  this  method  will  not  only  give  reliable  data  for  this 
range  of  values,  but  supply  the  information  in  near-real  time.  The 
precision  and  accuracy  of  the  data  are  being  evaluated. 

SUMMARY  AND  CONCLUSIONS 

Snow  cover  analysis 

Preliminary  analysis  of  the  Landsat  digital  data  using  the  GISS 
computer  algorithms  for  the  11  February  1973  scene  showed  that  the 
radiance  of  the  snow  cover/vegetation  varied  from  approximately  20 
mW/cm  sr  in  non-vegetated  areas  to  less  than  U mW/cm^^  sr  for  densely 
covered  forested  areas.  Comparison  of  the  digital  data  from  three  snow 
course  sites  in  the  Dickey-Lincoln  School  Lakes  area  with  the  radisince 
value  of  the  snowpeujk  at  these  sites  indicated  that  the  total  radiance 
of  th|  pixels  contained  in  the  snow  courses  varied  from  5.3U  tOp7.7U 
mW/cm  sr  with  the  average  radiance  value  being  6.U  + 0.6  mW/cm^^  sr. 

The  water  equivalent  of  the  snowpack  for  this  range  of  radiance  values 
was  approximately  2U.1  cm  (9*5  in.)  of  water.  Since  data  from  only 
three  snow  cowses  were  available,  the  correlation  between  radiance 
values  and  water  equivalent  of  the  snowpack  still  needs  to  be  tested. 
However,  if  the  relationship  hold  with  more  extensive  ground  truth  data 
in  the  Dickey-Lincoln  area,  it  is  anticipated  that  extrapolation  of 
radiance  values  for  snow  cover/vegetation  in  large  areas  of  the  watershed 
may  prove  useful  in  mapping  the  extent  of  snow  cover  vegetation  related 
to  a water  equivalent  value.  This  relationship  would  need  to  be  verified 
to  test  the  accuracy  of  using  satellite  imsigery  in  predicting  spring  run- 
off. 
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Table  VIII.  Typic®^-  data  from  the  tensiometer/transducer  interface  system. 

USA  CRREL  Test  Set 
Hanover,  New  Hampshire 
DCP  7110 

Calculations 


Date 

Time 

Tensiometer 
(cm  of  water) 

Error 

Water 

(volume) 

% Voids  filled 
with  air 

2/22 

1752 

-10.7 

27 

UO.O 

17.0 

2/23 

622 

-10.7 

07 

Uo.o 

17.0 

2/23 

:.8oi 

-22,5 

07 

Uo.o 

17.0 

2/2I4 

628 

-22.5 

07 

Uo.o 

17.0 

2/2U 

180U 

-22.5 

86 

Uo.o 

17.0 

2/25 

633 

-28.3 

07 

39.5 

17.5 

2/25 

81U 

-22.5 

07 

Uo.o 

17.0 

2/26 

638 

-28.3 

07 

39.5 

17.5 

2/26 

820 

-28.3 

07 

39.5 

17.5 

2/27 

6UU 

-3U.2 

07 

38.5 

19.7 

2/27 

829 

' -3U.2 

07 

38.5 

19.7 

3/1 

1832 

-40.1 

07 

37.8 

21.1 

3/2 

702 

-UO.l 

07 

37.8 

21.1 

3/2 

8UU 

-UO.l 

07 

37.8 

21.1 

3/3 

1703 

-51.9 

07 

36,0 

2U.8 

3/3 

2030 

-57.8 

80 

35. U 

26.1 

3/‘» 

856 

-57.8 

07 

35. U 

26.1 

3/6 

1721 

-81.3 

07 

30.2 

37.0 

3/6 

1859 

-81.3 

80 

30.2 

37.0 

3/7 

911 

-87.2 

07 

29.2 

39.0 

3/8 

735 

-93.1 

07 

28.8 

39.9 

3/8 

92U 

-99.0 

07 

27.9 

U1.8 

3/9 

7»*1 

-110.7 

07 

26.6 

UU.5 

3/9 

92U 

-110.7 

07 

26.6 

UU.5 

The  different  chemges  in  radiance  for  ^he  23  July  1973  ,'3  mW/cm 
sr)  arid  11  February  1973  (5-3^*  - 7.7**  mW/cm  sr)  Landsat  scenes  for  the 
three  snow  course  sites  were  attributed  to  the  presence  of  snow  and 
vegetation  changes  in  the  11  February  scene.  Multispectral  training 
signatures  for  four  snow  cOver/vegetation  classes  were  derived  from  the 
four- band  energy  values  for  the  11  Febrruary  1973  Landsat  scene.  These 
signatures  were  applied  to  the  digital  data  and  four  snow  cover/’vegeta- 
tion  classes  were  mapped. 


Landsat  digital  data  were  used  in  a wetlands  analysis  of  the 
Merrimack  River  est\iary  for  the  date  of  6 July  1976.  A multispectral 
signature  was  developed  for  a wetleinds  category  from  two  training  sites 
located  on  the  north  side  of  the  Merrimack  River  estuary.  Wetlands  were 
mapped  with  an  acctiracy  of  75^  when  compared  to  ground  truth  information. 

Mapping  of  flooded  areas 


Landsat  digital  data  were ‘also  used  to  delineate  flood  waters  in 
the  Franklin  Falls  reservoir  area.  New  Hampshire,  for  the  6 July  1973 
scene.  Low-water  reservoir  and  flood  water  stages  were  mapped  from 
grayscale  printouts  of  MSS  band  7 for  27  October  1972  and  6 July  1973, 
respectively.  Comparison  with  ground  truth  information  indicated  very 
good  agreement  with  the  Landsat  digital  data.  Approximately  60%  more 
water  was  observed  on  the  6 July  1973  scene  than  on  27  October  1972, 
which  was  as  expected. 


DCP  sensor  interface  development 


Snow  pillow  transducer  systems  for  measuring  the  water  equivalent 
of  the  snowpack  in  northern  Maine  were  interfaced  and  field  tested. 
Problems  with  temperature  sensitivity  of  the  transducers  were  encountered 
during  the  first  testing  period  (1975-76).  Problems  were  also  encountered 
diiring  the  second  testing  period  (1976-77 ) with  inadvertent  breakage  of 
the  transducers  and  questionable  data  from  a DCP  transmitter. 

A water  quality  monitor  interfaced  to  the  DCS  was  also  field  tested 
in  northern  Maine.  The  water  quality  data  compared  favorably  to  on-site 
chemical  analysis.  The  water  quality  data  obtained  will  be  used  as  base- 
line information  for  the  Dickey-Lincoln  School  Lakes  Project. 

A thermocouple  system  was  successfully  interfaced  and  field  tested  at 
Sugarloaf  Mountain,  Maine.  Temperature  data  from  the  sxarface  to  a depth 
of  30  m (lOO  ft)  were  transmitted  through  the  Landsat  DCS.  The  emplace- 
ment and  interfacing  techniques  developed  during  this  experiment  will  be 
used  for  future  installation  of  thermocouples  at  Alaskan  sites.  In 
addition,  this  temperature  measurement  system  could  be  used  in  reservoirs 
to  monitor  water  temperature  on  a daily  basis. 
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A tensiometer/transducer  system  was  successfully  interfaced  to  the 
Lands  at  DCS.  Laboratory  results  indicated  that  tension  as  low  as  10  cm 
(U  in.)  and  sis  high  as  900  cm  (35^  in.)  of  water  can  be  obtained. 
Presently  the  system  is  bfelng  tested  under  field  conditions. 

RECOMMENDATIONS 

The  following  recommendations  are  made  based  on  the  resitlts  of  the 
Landsat-2  investigation. 

1.  The  Landsat  CCT's  for  specific  su-eas  of  Interest  in  New  BngT 
should  be  furnished  to  the  user  on  a timely  basis  (within  three  days) 
for  use  in  evaluating  the  snow  cover  of  a watershed  during  critical 
spring  runoff  periods  (March  through  May ) . The  9~day  coverage  provided 
by  the  Landsat  satellites  is  not  adequate  for  the  operational  needs  of 
the  NED  Control  Center.  However,  in  times  of  critical  spring  flooding 
the  timely  use  of  the  available  Landsat  imagery  and  CCT's  would  be 
extremely  valuable  for  a regional  evaluation  of  the  extent  of  snow 
cover. 

2.  The  MSS  band  7 digltsLl  data  should  be  used  for  mapping  the 
extent  of  flood  waters  since  this  multispectrea  band  provided  the  :nost 
information  in  the  delineation  of  the  water/land  Ixjundary. 

3.  When  a mapping  accuracy  of  75/5  is  required,  the  Landsat  digital 
data  should  be  used  in  regional  evaluation  of  wetlands. 

remote  locations  the  Landsat  DCS  can  be  used  on  an  opera- 
tional basis  for  monitoring  hydrologic  and  environmental  parameters. 

FUTURE  PLANS 

The  snow  cover  analysis  work  will  continue  with  the  FTTS  Civil 
Works  sponsored  worx  unit  entitled,  "Snow  Cover  Analysis  in  New  England 
Using  Landsat  Digital  Data."  Site  selection  based  on  vegetation,  slope, 
aspect  and  elevation  will  be  accomplished  in  the  Dickey-Lincoln  School 
Lakes  project  area  and  other  selected  watersheds.  Ground  truth  measure- 
ments of  snow  depth  and  water  equivalent  at  the  selected  sites  will  be 
baken  in  conjunction  with  times  of  the  Lcuidsat  imagery  acquisition.  A 
meteorological  station  will  be  installed  to  obtain  data  on  local  climatic 
conditions.  As  available,  cloud-free  Landsat  CCT’s  will  be  acquired  for 
the  1977-78  winter  season  and  analyzed. 

Cloud-free  Landsat  CCT's  may  also  be  obtained  over  the  Sleepers 
River  Watershed  in  Danville,  Vermont,  for  the  1972-77  winter  seasons. 
Detailed  measxu'ements  of  the  snow  cover  are  available  at  this  site  from 
December  1968  to  present.  This  watershed  was  chosen  because  It  is 
hydrologically  representative  of  most  of  the  glaciated  upland  areas  in 
New  England  and  is  extensively  instrumented. 


Results  from  these  two  areas  will  be  developed  to  further  evaluate 
the  method  of  obtaining  the  water  equivalent  of  the  snow  cover  using 
Landsat  digital  data. 
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APPENDIX  B.  Computer  programs  used  to  convert  the  DCP  hexadecimal 
data  cards  into  real  numbers . 
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Figure  Bl.  Computer  program  to  decode  the  hexadecimal  data  cards  for 

the  snow  pillow  interface  located  at  Allagash  Falls,  Maine. 
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APPENDIX  C 


CIRCUIT  DESCRIPTION  OF  THERMOCOUPLE  INTERFACE 
Thermocouple  Interface 

The  thermocouple  vmit  was  designed  to  accept  inputs  from  up  to  112 
copper-constantan  thermocouples  which  are  arranged  in  l6  "banks"  of  7 
inputs.  Each  bank  of  inputs  is  recorded  with  an  individual  update  of 
the  data  collection  platform  (DCP).  There  were  28  thermocouples  (four 
banks)  used  at  the  Sugeu-loaf  Moxmtain  installation. 

The  reference  thermocouple  Junction  is  compensated  to  within  +0.3®C 
over  a temperature  range  of  -50°  to  +Uo°C  using  a combination  of  five 
thermistors.  The  normal  measurement  temperature  range  is  -3^°  to  +32°C 
with  a resolution  of  +0.25°C  or  10  pV.  It  is  possible  to  trade  range 
for  resolution,  or  vice  versa,  by  selecting  different  groups  of  data 
bits  from  the  analog-to-digiteLl  converter. 

The  thermocouple  unit  is  designed  to  operate  with  a LaBarge  Elec- 
tronics Convertible  Data  Collection  Platform  (C/DCP)  which  has  a memory 
capability.  Power  requirements  are  +12  V (nominal)  at  0.5  ampere  and  5 V 
at  1.0  milliampere.  The  12-V  power  is  applied  only  while  the  C/DCP  is 
acquiring  new  data. 

Circuit  Description 

The  electronic  components  are  arranged  on  nine  circuit  cards  that 
are  identified  as  follows  (Fig.  Cl): 

Input  multiplexer,  group  I (MX  I) 

Input  multiplexer,  group  II  (MX  II) 

Input  multiplexer,  group  HI  (MX  III) 

Input  multiplexer,  group  IV  (MX  IV) 

Amplifier  and  channel  multiplexer  (AMPL/CHANNEL  MUX) 
Analog-to-digital  converter  (A/D  CONV) 

Power  supply  and  reference  Junction  compensator  (P.S./REF  JCT  COMP) 
Latch  (LATCH) 

Interface  (INTERFACE) 

Thermocouple  signals  are  routed  through  the  input  multiplexer,  con- 
sisting of  cards  MX  I,  MX  II,  MX  III  and  MX  IV  (Fig.  C2,  C3,  CU,  C5). 
Seven  input  signals  are  read  diu'ing  one  C/DCP  update  sequence.  The  l6 
banks  are  numbered  0 through  15.  When  an  update  occurs,  the  next  high- 
est numbered  bank  is  read.  Thus,  after  l6  updates  a total  of  112  inputs 
have  been  sampled. 
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Each  of  the  seven  inputs  of  the  selected  bank  occupies  an  8-bit 
word  of  C/DCP  data  and  the  inputs  are  desifcnated  as  channels  1 through 
7.  The  remaining  word,  corresponding  to  channel  0,  is  the  number  of  the 
bank  selected.  Four  of  the  eight  bits  of  this  word  are  required  for 
this  process  suid  the  remaining  four  are  always  zero.  There  is  an  avail- 
able analog  input  for  channel  0.  It  is  used  only  during  bench  testing 
and  adjusting  and  has  no  corresponding  C/DCP  data  word. 

The  amplifier  and  chsuinel  miiltiplexer , consisting  of  one  card 
(.\MPL/ CHANNEL  MUX),  scans  the  seven  charmels  and  amplifies  their  signals 
to  a suitable  level  for  digitizing  (Fig.  C6).  This  module  also  controls 
the  associated  analog-to-digited.  converter  (A/D  CONV)  and  provides 
signals  (called  latch  strobes)  which  control  the  treinsfer  of  the  digi- 
tized data  to  the  latch  card  (Fig.  Cj). 

The  A/D  CONV  is  a prepackaged  unit  mounted  on  a separate  card  (Fig. 
C7).  Although  it  has  a 12-bit-plus-sign  capability,  only  7 bits  and  the 
sign  appear  in  the  DCP  data.  The  A/D  CONV  has  a differential  input,  in 
addition  to  the  normal  analog  signal  input,  and  receives  the  output  of 
the  reference  Junction  compensator. 

The  system  power  supply  and  reference  Junction  compensator  modules 
(P.S./REF  JCT  COMP)  are  mounted  on  one  card  (Fig.  C8).  The  reference 
Junction  compensator  generates  a voltage  which  matches  that  produced  by 
the  reference  Junction  itself  over  the  temperature  range  of  -50°  to 
+40°C.  The  compensating  voltage  is  introduced  to  the  A/D  CONV  at  the 
same  level  as  the  amplifier  output.  The  compensating  voltage  is  the 
reference  Junction  voltage  multiplied  by  the  gain  of  the  amplifier. 

The  latch  card  (LATCH)  contains  eight  8-bit  latches  to  accept 
parallel  data  (Fig.  C9).  A strobe  signal  is  required  to  enter  data  and 
an  enable  signal  is  needed  to  make  data  available  at  the  output.  If  the 
appropriate  signal  is  not  present,  the  device  appears  to  be  an  open 
circuit.  Thus  the  input  and  output  ports  can  operate  on  common  busses. 
The  8-bit  parallel  data  from  the  A/D  CONV  is  bussed  to  all  of  the  latch 
inputs.  A strobe  signal  from  the  channel  multiplexer  operates  the  input 
of  each  latch  in  succession,  according  to  the  analog  channel  being 
sampled. 

The  data  input  to  the  C/DCP  is  by  way  of  a l6-bit  parallel  bus. 

The  8-bit  latches  are  paired,  so  that  the  output  lines  of  latch  pairs 
(0-1,  2-3,  l*-5  and  6-7)  each  form  a l6-bit  data  word. 

The  enable  signals  for  all  latches  (except  0 and  l)  are  received 
from  the  C/DCP.  The  enable  signals  are  the  triggering  signals  that 
would  normally  be  used  to  start  the  punch  motors  in  the  analog-to- 
digital  recorders  such  as  are  operated  by  the  U.S.  Geological  Survey. 
These  are  labeled  "Punch  motor  no . 1 , " "Punch  motor  no . 2 , " etc . , in  ed.1 
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Figure  C9.  Latch  (LATCH). 


instructions  and  diagrams  (Fig.  Cl).  Punch  motor  no.  1 initiates  a scan 
hy  the  channel  multiplexer.  Punch  motor  signals  2,  3 and  k enable  latch 
pairs  2-3,  ^-5  and  6-7,  respectively,  to  be  started.  Since  Punch  motor 
no.  1 arrives  before  latches  0 and  1 can  be  supplied  with  data,  the 
input  strobe  signal  for  latch  7 also  serves  as  the  output  enable  signal 
for  latches  0 and  1. 


Latches  1 through  7 carry  thermocouple  data,  which  consist  of  7 
magnitude  bits  and  one  sign  bit.  Latch  0 carries  a U-bit  word  to  indi- 
cate which  of  the  l6  banks  of  inputs  is  being  scanned.  The  remaining  U 
bits  in  latch  0 are  not  presently  utilized. 

The  l6-bit  output  from  the  latches  goes  throxigh  a transistor  switch 
array  (INTERFACE)  to  the  C/DCP  (Fig.  CIO).  The  transistors  invert  the 
polarity  of  the  data  signals  and  also  simulate  the  grounding-type  signal 
that  would  come  from  an  analog-to  digital  recorder. 
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Figure  CIO.  Interface  (INTERFACE). 
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The  primary  power  supply  is  12  V DC,  supplied  from  the  battery  that 
powers  the  C/DCP,  and  is  turned  on  dxiring  an  update  sequence  by  a 
semiconductor  switch  in  the  C/DCP.  The  switched  12-V  supply  powers  a 
+15-V  and  a +5-V  supply  which  sure  electrically  isolated  from  each  other 
to  avoid  current  loops.  The  bank  select  counter  and  certain  parts  of 
the  control  portion  of  the  amplifier/channel  multiplexer  require  +5  V 
continuously.  The  voltage  comps  from  a regulated  supply  within  the 
C/DCP. 

Operating  Sequence 

The  operating  sequence  is  as  follows : 

1.  The  C/DCP  begins  an  update  sequence,  by  applying  a 12-V  input  to  the 
^15-V  and  +5-V  power  supplies,  and  initializing  a "Punch  motor  no.  l" 
trigger  pulse  to  the  control  portion  of  the  amplifier  and  cheuinel  multi- 
plexer. This  resets  the  channel  multiplexer  to  channel  0.  At  the  same 
time,  the  channel  multiplexer  sends  a clock  pulse  to  the  bank  select 
counter,  which  is  on  MX  I,  advancing  it  by  the  count  of  one. 

2.  The  channel  multiplexer  scans  the  seven  channels  of  the  bank 
selected.  It  remains  on  each  channel  for  about  0.5  second,  diuring  which 
time  a pulse  is  sent  to  the  A/D  CONV  to  start  a conversion.  When  the 
conversion  is  complete,  a pulse  (latch  strobe)  is  sent  to  the  proper 
latch  to  enter  the  data.  This  process  is  repeated  for  all  channels.  On 
channel  0,  a conversion  is  performed,  but  channel  0 data  are  not  entered 
into  a latch.  Instead,  the  state  of  the  beuik  select  counter  is  entered. 
Thus,  upon  completion  of  a scan,  latch  0 contains  the  bank  number,  while 
latches  l through  7 contain  temperat\ire  data. 

3.  The  data  in  latches  0 through  8 are  loaded  into  the  C/DCP  memory. 

As  previously  described,  the  C/DCP  "P\mch  motor"  signals,  together  with 
the  strobe  signal  for  latch  7,  enable  the  latches,  two  at  a time,  to 
operate. 

U.  Upon  completion  of  the  update  sequence,  which  requires  about  90 
seconds,  the  C/DCP  removes  the  12-V  power  from  the  lanit.  The  C/DCP 
continues  to  supply  5-V  power  to  the  bank  select  counter  and  certain 
other  circuits  which  must  remain  energized. 


In  accordance  with  letter  from  DAEN-RDC,  DAEH-ASI  dated 
22  July  1977 » Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a facsimile  catalog 
card  in  Library  of  Congress  MABC  format  is  reproduced 
below. 
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